Electronic structure of BAs and boride III— V alloys 
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Boron arsenide, the typically-ignored member of the III-V arsenide series BAs-AlAs-GaAs- 
InAs is found to resemble silicon electronically: its F conduction band minimum is p-like (ris), 
not s-like (ri c ), it has an Ai c -like indirect band gap, and its bond charge is distributed almost 
equally on the two atoms in the unit cell, exhibiting nearly perfect covalency. The reasons for 
these are tracked down to the anomalously low atomic p orbital energy in the boron and to 
the unusually strong s-s repulsion in BAs relative to most other III-V compounds. We find 
unexpected valence band offsets of BAs with respect to GaAs and AlAs. The valence band 
maximum (VBM) of BAs is significantly higher than that of AlAs, despite the much smaller 
bond length of BAs, and the VBM of GaAs is only slightly higher than in BAs. These effects 
result from the unusually strong mixing of the cation and anion states at the VBM. For the 
BAs-GaAs alloys, we find (i) a relatively small (~3.5 eV) and composition-independent band 
gap bowing. This means that while addition of small amounts of nitrogen to GaAs lowers the 
gap, addition of small amounts of boron to GaAs raises the gap (ii) boron "semi-localized" 
states in the conduction band (similar to those in GaN-GaAs alloys), and (iii) bulk mixing 
enthalpies which are smaller than in GaN-GaAs alloys. The unique features of boride III-V 
alloys offer new opportunities in band gap engineering. 



I. INTRODUCTION 

With the advent of state-of-the-art techniques for 
growing semiconductor alloys on common substrates such 
as GaAs, silicon, and germanium, semiconductor com- 
pounds which previously were very difficult to synthesize 
are now routinely achieved. Techniques such as metal- 
organic chemical vapor deposition (MOCVD), molecular 
beam epitaxy (MBE), and pulsed laser ablation (PLD) 
have provided the opportunity to synthesize and study 
a large number of nitride, phosphide, and antimonidc 
semiconductor alloys. Of particular recent interest are 
alloys of a wide gap semiconductor (e.g. nitrides) with 
a "conventional" III-V semiconductor because of their 
promise in optical applications. Two diverging scenarios 
were considered: (i) using a significant amount (10%- 
30%) of the wide gap component to shift the alloy band 
gap to the blue (e.g., adding ^20% GaN to InN) for light- 
emitting diode or laser applications, and (ii) using a small 
amount of the wide gap semiconductor to shift the alloy 
band gap to the red (e.g., adding 1-3% GaN to GaAs) for 
photovoltaic applications. The latter effect occurs natu- 
rally if the band gap bowing parameter b is larger than 
the difference of the band gaps of the constituents (e.g., 
ZnS-ZnTe; GaAs-GaN). In this case addition of small 
amounts of the wide gap components acts to initially 
lower the band gap of the small gap component. For ex- 
ample, one can achieve the technologically desired 1 eV 
gap if one adds nitrogen to GaAs or to InGaAs. 

When boron is substituted into GaAs, it can go to ei- 
ther a gallium site or an arsenic site. Normally boron 



prefers isovalent substitution on the gallium sitejdotl 
which is the case we study here. In the other case when 
boron goes to the arsenic site (a boron "antisite" defect) , 
the boron acts as an acceptor and this -antisite defect 
has been the subject of numerous studies. qEtdQeI Growth 
conditions determine whether boron goes to the gallium 
site as an isovalent substition or to the arsenic site as an 
acceptor. For example, B,4 S antisite defects are known 
to occur in Ga-rich samples of GaAs grown by the liquid 
encapsulated Czochralski (LEC) method, but for GaAs 
crystals taken from stoichiometric or As-rich melts, elec=. 
trically active boron or boron complexes are not found, tl 
indicating that the boron atoms have substituted isova- 
lently to the gallium sites. Not as much is known about 
the epitaxial growth conditions under which isovalent or 
antisite boron incorporation occurs. It is reasonable to 
suppose that B^ s antisite defects will be more likely un- 
der Ga-rich conditions and that As-rich growth condi- 
tions will lead to isovalent boron incorporation, similar 
to the case for LEC-grown GaAs. In this study, we focus 
on isovalent BGaAs alloys where boron occupies gallium 
sites. 

In this paper we will explore BAs as an alternative 
to GaN as a wide gap partner for alloying conventional 
III-Vs such as AlAs and GaAs. Very little is known 
about BAs: As early as 1966, Ku mentioned the possi- 
bilities of BAs-GaAs as aJapside III-V alloy with po- 
tentially useful properties .1190 However, difficulties in 
fabricating. BAs <ff pimple solid solutions of zinc-blende 
BAs-GaAsBBEayy prior to the development of cur- 
rent epitaxial techniques have prevented a careful study 



of the alloy properties. Now that epitaxial techiques have 
eased these difficulties somewhat, boridenSeniicQaductor 
alloys are generating renewed interest BtallSIIZHIa There 
are a number of theor etica l studios that .examine differ- 
ent aspects of pure BAstBE3E2lE3'E3ESE3E3 but theoretical 
studies of the boron alloys are lacking. 

Important questions that one would want to answer 
include: 

(i) How will alloying with BAs affect the band gaps 
and other properties of GaAs or AlAs? 

(ii) To what extent does BAs fit the well-known trends 
in band gaps, band offsets and bonding patterns of the 
arsenide III-V series BAs-AlAs-GaAs-InAs? 

(iii) Will the band gap bowing of B^Gai-^As alloys 
be as anomalous as the (very large and composition- 
dependent) bowing in GaAsi-^N^? 

(iv) Will cation substitution by boron lead to unusual 
wavefunction localization effects found to exist for anion 
substitition by nitrogen?E3c3'E£l 

Two features of boron make boride semiconductors 
fundamentally different from common III-V or II-VI 
semiconductors. The first is that, like nitrogen, boron 
is in the first row of the periodic table and has deep p 
orbitals and a small atomic size. The second feature is, 
unlike nitrogen, boron has a low electronegativity. This 
leads to highly covalent compounds, unlike nitride semi- 
conductors, which have a strong ionic character. This 
paper examines (i) zinc-blend BAs and its place in the 
III- As family of semiconductors and (ii) boron substitu- 
tion of gallium in GaAs, including alloy bowing, band 
offsets and mixing enthalpies. Our main findings are: 

Zinc-blende BAs: Surprisingly, we find that, electroni- 
cally, BAs resembles silicon rather than other III-V semi- 
conductors. Similar to silicon and in contrast to most 
III-Vs, the lowest Brillouin-zone center conduction band 
of BAs has p symmetry (ri5 C ) rather than s symmetry 
(Tic), and, like silicon, its total valence charge density 
shows almost symmetric distribution of charge around 
the two atoms in the unit cell. The reasons for the silicon- 
like conduction band ordering in BAs are: (i) the small 
repulsion of the bonding and antibonding p states due 
to the low p orbital energy of boron, as well as the un- 
usual hybridization of both cation and anion p states at 
the VBM, and (ii) the repulsion of the cation and anion 
s states which is much stronger in BAs than in AlAs, 
GaAs, and InAs. As a result of the p-p hybridization 
(covalent bonding) mentioned in (i), we also find that 
the valence band offset of BAs relative to other members 
of the III- As family is unusually high. 

BAs-GaAs alloys: (i) The band gap bowing is rela- 
tively small (~3.5 eV) and composition- independent, in 
stark contrast to GaN-GaAs alloys. Because of this small 
bowing, the addition of BAs to GaAs increases the gap, 
thus, unlike nitrogen, addition of boron into GaAs or In- 
GaAs will not lead to the desired 1 eV material, (ii) The 
lower energy conduction band states are "semi-localized" 
states around the boron atoms, e.g., the CBM is strongly 
localized near the boron but extended at longer distances 



while the VBM is completely delocalized. (iii) The bulk 
mixing enthalpy of BAs in GaAs is much lower than that 
of GaN in GaAs, indicating that the bulk solubility of 
boron in III-V compounds may be higher than that of 
nitrogen and thus higher composition ranges may be pos- 
sible with the boride alloys. These findings indicate that 
boride III-V alloys provide new opportunities in band 
gap engineering. 



II. METHODS OF CALCULATION 
A The LAPW calculations 

We used density functional— theory within the local 
density approximation (LDA),E3 as implemented by the 
full-potential linearized-augmented-plajiewaves (LAPW) 
methodtilcS (wien97 implementationEa) . The exchange- 
correlation potential of Perdew and Wang was used.E-3 In 
the calculations with less than 32 atoms, the plane wave 
kinetic energy cutoff for the expansion in the intersti- 
tial region was 16 Ryd (approximately 130 basis function 
per atom). The muffin tin (MT) radii were 1.65 bohr 
for boron and 2.2 bohr for arsenic, aluminum, gallium, 
and indium. In the large supercell calculations (32 or 
more atoms), a slightly smaller plane wave kinetic en- 
ergy cutoff 13 Ryd (approximately 120 basis functions 
per atom) and a larger boron MT radius of 1.8 bohr (2.1 
for arsenic, aluminum, gallium, and indium) was used to 
ease the computational burden of the larger cells. Our 
convergence studies indicate that the error in the indi- 
vidual eigenstates is less than 5 meV for the valence 
and lower conduction bands. The calculations were run 
until the variation in the total energy between several 
self-consistency cycles was < 10~ 5 Ryd. The experimen- 
tal lattice constants were used in all the calculations of 
the individual compounds. The experimental lattice con- 
stants are 4.777 A, 5.660 A, 5.653, and 6.058 A for BAs, 
AlAs, GaAs, and InAs, respectively. 

The k point mesh used in the calculation of the simple 
binary compounds (zinc-blende BAs, AlAs, GaAs, and 
InAs) was a 4 x 4 x 4 mesh of Monkhorst and Pack spe- 
cial points (10 points in the irreducible wedge of the Bril- 
louin zone [BZ] ) Ei The superlattice calculations for the 
valence band offsets as well as the supercell calculations 
for the alloy studies used k point meshes equivalentEa to 
the 4x4x4 mesh used in the calculation of the sim- 
ple zinc-blende binary compounds. Using equivalent k 
point meshes is particularly important for calculations 
such as enthalpies of formation in order to eliminate un- 
certainties due to the statistical errors of different k point 
meshes. Thus, one can use much smaller k point meshes 
to achieve the required accuracies than would otherwise 
be necessary. 



B Partial DOS, band characters, and valence 
charge density 

It is useful to analyze the orbital character of differ- 
ent states. The band character (or orbital population) 

Q" (e, k) is the £-th angular momentum component of 
the charge due to wavefunction ip(e, k) enclosed in a 

sphere £Imt of radius Rm T about atom a 



Q ( e a) (e,k) = 



|P^(e,k,r)| 2 dr 



(1) 



where Pi is an angular momemtum projection operator 
with origin at site a. Because the interstitial region out- 
side the muffin tin spheres is excluded when the band 
characters are calculated, the muffin tin radii were chosen 
to match the rationalized tetrahedral radii of Phillips.LJ 
The partial DOS are determined by integrating 
Q e a (e,k) over all k in the Brillouin zone (BZ); 



<\e) 



Q< Q) (e,k)dk. 



(2) 



BZ 



The angular-momentum decomposition of the total va- 
lence charge can then be determined by summing Ni(e) 
over the valence bands; 



VB„ 



vb„ 



Nt(e)de. 



(3) 



Valence charge densities: The charge density is con- 
structed from the highest Nb occupied bands as 



Pval(r) 



N B 

E 



|V> n ,k(r)| 2 dk. 



(4) 



BZ 



For plotting p va \ we use Nb = 4. Thus, cation d bands 
are not included in the construction. 

The "valence deformation density" Ap va i(r) describes 
the difference between the solid-state density and a su- 
perposition of densities of spherical atoms in their ground 
state, 

Ap val (r) = p va i(r) - p S up(r)- (5) 

The valence charge difference is defined as the difference 
between the crystal valence density (including the cation 
d bands, when occupied) and a superposition of free, 
spherical LDA atomic densities arranged in the config- 
uration of the crystal. The atomic configuration of the 
free atoms are d 10 s 2 p 
anions. 



for the cations and s 2 p 3 for the 



C Band offsets 



The offset AE V (AX/BY) between the valence band 
maxima of two semiconductor compounds AX and BY 



TABLE I: Convergence test for the valence band offsets 
of BAs/GaAs, BAs/AlAs, and AlAs/GaAs using two dif- 
ferent superlattice periods. Units are eV. Band offsets 
were calculated using both anion and cation core states 
[see Eq. (0)]. Note that while for AlAs/GaAs a similar 
offset is obtained when either anion or cation core levels 
are used, for the borides only anion core levels provide a 
rapidly convergent band offset with respect to the super- 
lattice period n. 



Superlattice 


Using 
f s anion 


Using 
2s anion 


Using 
Is cation 


(BAs) 2 /(GaAs) 2 
(BAs) 4 /(GaAs) 4 


0.18 
0.19 


0.11 
0.12 


0.83 
0.58 


(BAs) 2 /(AlAs) 2 
(BAs) 4 /(AlAs) 4 


-0.38 
-0.39 


-0.45 
-0.47 


0.58 
0.17 


(AlAs) 2 /(GaAs) 2 
(AlAs) 4 /(GaAs) 4 


0.50 
0.51 


0.50 
0.51 


0.41 
0.47 



forming a heterostructure is calculated using a method 
similar to that used in photoemission spectroscopy!^! 



E V (AX/BY) = AE^, 



AE^ + AE^ BY , (6) 



where 






E BY 



E$ x 

Eq, 



(7) 

are the energy seperations between the core levels (C and 
C) and the valence band maximums for the pure AX and 
BY compounds. The third term in equation (o), 

AE A C X C IBY = E*Y - E^ (8) 

is the difference in the core level binding energy be- 
tween the two compounds AX and BY in the AX/ BY 
heterojunction. To calculate the "natural band offset" , 
A£^£ and A£^£, are calculated for AX and BY at 
their cubic, equilibrium lattice constants. The core level 
difference AE C c , is obtained from calculations of 
(AX) n /(BY) n superlattices in the (001) direction. The 
period n needs to be large enough so that AA-like and 
By-like properties can be identified far from the inter- 
face. For the core states (denoted C and C"), we use 
the anion Is states. Table [fl shows the dependence of the 
valence band offset calculations on the period of the su- 
perlattice and on which core states are used. Because of 
the large size mismatch of BAs and GaAs, at least n = 4 
is required to sufficiently converge AE C c , . We also 
find that while a similar offset is obtained for AlAs/GaAs 
when either anion or cation core levels are used, for the 
borides only anion core levels provide a rapidly conver- 
gent band offset with respect to the superlattice period 
n. 



TABLE II: Definition of the supercells used in this study and equivalent k points. Lattice vectors are given in units 
of do and equivalent k points are given as fractions of the reciprocal lattice vectors. 



System 



AC, BC zinc-blcndc 



Lattice vectors 



(1/2, 1/2, 0) 
(1/2, 0, 1/2) 
(0, 1/2, 1/2) 



Equivalent k points 



(0, 0, 1/8) 
(0, 0, 3/8) 

(0, 1/8, 3/4) 

(0, 1/8, 1/4) 

(0, 1/8, 1/2) 

(0, 1/4, 5/8) 

(1/8, 1/4, 1/2) 

(0, 1/4, 3/8) 

(0, 3/8, 1/2) 

(1/8, 3/8, 5/8) 

(0, 1/8, 1/8) 

(0, 1/8, 3/8) 

(0, 3/8, 1/8) 

(0, 3/8, 3/8) 

(1/8, 1/4, 1/8) 

(1/8, 1/4, 3/8) 

(1/8, 1/2, 1/8) 

(1/8, 1/2, 3/8) 

(1/4, 3/8, 1/8) 

(1/4, 3/8, 3/8) 

(3/8, 1/2, 1/8) 

(3/8, 1/2, 3/8) 



Relative weight 



(iC)i/(BC)i (100) supcrlattice 



(1/2, 1/2, 0) 

(1/2, -1/2, 0) 

(0, 0, 1) 



(AC) 2 /{BC) 2 (100) supcrlattice 



(1/2, 1/2, 0) 

(1/2, -1/2, 0) 

(0, 0, 2) 



(0, 1/8, 1/4) 
(0, 3/8, 1/4) 
(1/8, 1/4, 1/4) 
(1/8, 1/2, 1/4) 
(1/4, 3/8, 1/4) 
(3/8, 1/2, 1/4) 



(AC)i/{BC)i (100) supcrlattice 



(1/2, 1/2, 0) 

(1/2, -1/2, 0) 

(0. 0, 4) 



(0, 1/8, 1/2) 
(0, 3/8, 1/2) 
(1/8, 1/4, 1/2) 
(1/8, 1/2, 1/2) 
(1/4, 3/8, 1/2) 
(3/8, 1/2, 1/2) 



AiBiCs fee supcrccll 



(1, 1,0) 
(1,0, 1) 
(0, i, i) 



(0, 0, 1/4) 
(0, 1/4, 1/2) 



AiSigCig bec supcrccll 



irnr-i) 
(1,-1, 1) 
(-1, 1, 1) 



(1/8, 1/8, 1/8) 
(1/8, 5/8, 5/8) 



A1B31C32 simple cubic supcrccll 



AgBsCie SQS16 



(2, 0, 0) 

(0, 2, 0) 
(0, 0, 2) 



(1,-1,2) 

(1,-1-2) 

(1/2, 1/2, 0) 



(1/4, 1/4, 1/4) 



(0, 1/2, 0) 

(1/2, 0, 0) 

(1/4, 3/4, 1/8) 

(0, 1/2, 1/4) 

(1/2, 0, 1/4) 

(1/4, 3/4, 3/8) 

(0, 1/2, 1/2) 

(1/2, 0, 1/2) 



D Alloy enthalpies of mixing, alloy bowing, and 
equivalent k points 

The enthalpy of mixing for an alloy A\- X B X of two 
components A and B is the difference in energy between 
the alloy and the weighted sum of the constituents; 



AH{x) = E ApBq - [(1 - x)E A + xE B ] 



(9) 



where x = p/(p + q). In the mixing enthalpy calculations, 
the alloy lattice constant a{x) is taken as the linear aver- 
age of the experimental values of the constituents. The 
experimental lattice constants of BAs and GaAs are 4.777 
A and 5.653 A, respectively, whereas the calculated lat- 
tice constants are 4.740 A and 5.615A. Any free internal 
coordinates in the alloy structure were relaxed using the 



quantum mechanical forces so that residual forces on the 
ions were less than 1 mRyd/bohr. 

The band gap bowing parameter b defined by 



E g (x) = E g (x) - bx(x - 1), 



(10) 



where E g (x) is the weighted linear average of the indi- 
vidual band gaps of A and B. The bowing parameter b 
represents the deviation of the band gap from this av- 
erage. Note that although there are LDA errors in the 
calculated gaps of the A X B\_ X alloy as well as the pure 
constituents A and B, to lowest order, the LDA error in 
the bowing b cancels out. 

Considering Eqs. (Q) and (|10|), we see that one needs to 
converge the k representation for a compound A p B q as 
well as for the elemental constituents A and B. The stan- 



TABLE III: All possible pair configurations in the 64 
atom cell. 



Shell number 
(nth fee neighbor) 



Lattice vector Number of 

[ti,V, w] (unit ao/2) equivalent pairs 



[0, 1, 1] 


12 


[2, 0, 0] 


3 


[2, 1, 1] 


12 


[2, 2, 0] 


3 


[2, 2, 2] 


1 



dard way of accomplishing this is to increase the num- 
ber of k points in all three systems until convergence is 
obtained. The disadvantage of this approach is that it 
requires absolute k point convergence for A and B, and 
separately for A p B q . A better approach-is to take ad- 
vantage of relative k point convergence!!!] The idea is 
to sample the Brillouin zone equivalently for A, B and 
A p B q . This could be done by considering A p A q , B p B q , 
and ApB q as isostructural solids and sampling the Bril- 
louin zone of each equally. Then, any relative fc-point 
sampling error cancels out. This is called the method of 
equivalent k points™ In practice, one does not have to 
calculate A p A q and B p B q but instead can calculate A 
and B in their primitive unit cells using suitably folded- 
in k points. Equivalent k points for the unit cells in this 
paper are given Table [n|. 



E Choice of supercells 



The calculations of AH [Eq. (|)] and b [Eq. @] 
require supercells. We use BiGarAss, BiGai5Asi 6 , 
BiGa3iAs32, and B2Ga3oAs32. The lattice vectors defin- 
ing the supercells are given in Table O. The SQS16 su- 
percell is a "special quasirandom structure" — a periodic 
structure with rather small unit cell whose lattice sites 
are occupied by A and B atoms so as to mimic the atom- 
atom correlation functions of rrmch larger A\^ X B X super- 
cells with random occupations SB 

In the calculations for the band gap bowing of 
B^Gai-^As alloys, 64 atom, simple-cubic unit cells were 
used for both the 3% and 6% boron alloys. In the case 
of the 3% alloy, there is one boron atom in the super- 
cell, but for the 6% alloy, there are two boron atoms in 
the supercell. For this case, the band gap was deter- 
mined by taking the weighted average of the gaps for 
the 5 symmetrically-inequivalent configurations (given in 
Table HI) of two boron atoms in the 64 atom supercell. 
These 5 pairs are the first through fourth neighbor pairs 
in an fee lattice, as well as the sixth neighbor. (In the 
64 atom cell, the fifth fee neighbor is equivalent to the 
first.) Using the same 64 atom unit cell for the alloys 
in the bowing calculations eliminated band gap differ- 
ences that can occur due to the k-point folding relations 
of different supercells, as discussed by Bellaiche et al. in 
Rcf. |2l. 



TABLE IV: Parameters used in the VFF model. d° ond is 
the equilibrium (unstretched) bond length, a is the bond 
stretch term, /? is the bond bending term, and a is the 
stretching-bending term. The units of a, j3, and a are 
10 3 dyne. 



System 


rfbond(A) 


a 


P 


a 


For band gap bowing 


BAs 

GaAs 




2.0685 
2.4480 


76.26 
32.15 


19.12 
9.37 


-9.12 
-4.10 


For mixing enthalpy 


BAs 
GaAs 

GaN 




2.0685 
2.4480 
1.9520 


76.26 
41.19 
96.30 


22.2 
8.94 
14.8 


— 



F Valence force field model 

In two cases, the total energies or the internal struc- 
tural parameters were determined using a classical 
force |S<|W| model. This generalized valence force field 
modelcffl (VFF) consists of bond stretching, bond bend- 
ing, and bending-stretching springs with force constants 
<7, /3, and a, respectively. The force constants and the 
equilibrium (no force) bond lengths are given in Table 
|yL The two cases where the VFF model was used are: 

(1) The VFF model was used to model the mixing en- 
thalpies of boron or nitrogen in G aAs in the regime of 
very low concentrations (Sec. |V C| ). Because well tested 
force constants exist for the "conventional" VFF model 
(a, (3; no a) for the case of GaAs/GaN, we also used 
a two-parameter (a = 0) model for BAs/GaAs for the 
sake of consistency with previous calculations. The pa- 
rameters for BAs/GaAs were tested by comparing the 
enthalpies of formation for structures at higher concen- 
trations calculated using first principles, as given in Ta- 
ble M 

(2) Full relaxation using first principles' forces of all 
64 atom supercells us ed in the band gap bowing calcula- 
tions reported ua-|Sec. VB is computionally prohibitive so 
the generalizedcll VFF model was used to determine in- 
ternal coordinates. In one case (BiGa3iAs32), the VFF- 
and LAPW-determined coordinates were compared and 
found to be very similar. The average difference between 
the two methods for the Ga-As bonds was < .005 a.u. 
(0.1%). The largest discrepency was near the boron site 
where the VFF predicts a B-As bond length that is 0.01 
a.u. (1.5%) shorter than the LAPW bond length. Using 
the VFF-determined internal coordinates of the super- 
cell would be acceptably accurate only if the band gap is 
very close to that when the first-principles' coordinates 
are used. In the case of the BGa3iAs32 supercell just 
mentioned, the band gap difference was only 10 meV, 
justifying the use of the VFF. 



TABLE V: Enthalpies of formation (LDA and VFF) of ordered (001) superlattices of BAs/GaAs and the mixing 
enthalpies of B^Gai-^As random alloys [Eq. (B)]. All the AH values are positive. 



System 


Boron concentration 


Ai7(meV/atom) 


A#(meV/atom) 






(LAPW) 


(VFF) 




Ordered superlattices: 






(BAs)i/(GaAs)i (001) 


50% 


186 


193 


(BAs) 2 /(GaAs) 2 (001) 


50% 


171 


192 


(BAs) 4 /(GaAs) 4 (001) 


50% 


161 


191 


(BAs) 00 /(GaAs) 00 (001) 


50% 


148 


— 


(BAs)i/(GaAs) 3 (001) 


25% 
Random alloys: 


132 


174 


BGa 7 As 8 (fee) 


12.5% 


73 


79 


BGaisAsie (bee) 


6.25% 


22 


30 


BGa3iAs32 (simple cubic) 


3.13% 


9.5 


16 



III. BAs AND THE III-As FAMILY 

A Expectations from atomic physics 

The LDA valence orbital energies for the first four el- 
ements of column III are plotted in Fig. fil The lower 
value of the atomic p states of B relative to the atomic 
p states for Al, Ga, and In (which are all quite similar) 
are notable. This lower value can be understood from 
the fact that, as a first row element, the 2p states in B 
need not be orthogonal to lower p states. In general, the 
cation s energy should increase going down the column 
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FIG. 1: Atomic s and p orbital energies in column III 
of the Periodic Table. Calculated fully rclativistically 
within the local density approximation. 



but there is a kink when moving from Al to Ga due to the 
introduction of imperfectly screened d states in Ga. The 
fact that the p states in B lie lower in energy, and hence 
closer to the As Ap states, will lead to a much stronger 
hybridization of p-like states in BAs relative to the other 
III-As systems. 

The consequence of this is shown in Fig. which shows 
the total valence charge, g^ ot of angular momentum £ [Eq. 
@], enclosed in spheres having the tetrahedral radii of 
Phillips .til It is clear that the occupation of the B and 
As p-like states in BAs is much higher than in the other 
systems (indicating that p-p hybridization is much more 
pronounced in BAs) and also that there is more charge 
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FIG. 2: Total valence charge [Eq. (0)] enclosed in spheres 
having the tetrahedral radii of PhillipsrH boron 0.853 
A, aluminum 1.230 A, gallium and arsenic 1.225 A, and 
indium 1.405 A. 



12.0 




FIG. 3: Relativistic energy band structure of BAs. The 
VBM is taken to be zero and the lattice constant is 4.777 

A. 



around the boron than the other cations. Thus, the dis- 
tinction between "cation" and "anion" is not so clear cut 
in BAs, where both atoms share similar charge. This is 
also reflected in the relative ionicities in the group. Ac- 
cording to the Phillips scalecHl the relative ionicities are 
.002, .274, .310, and .357 for BAs, AlAs, GaAs, and InAs, 
respectively. As a result of the relative ionicities, the 
bonding in BAs is almost completely covalent whereas 
the other three members of the group have a signifi- 
cant ionic component that increases in the order AlAs 
— ► GaAs — ► InAs. We will see next how these expecta- 
tions are reflected in the electronic structure of BAs. 



Energy bands, band characters, and densities of 
states in zinc-blende BAs 



The energy bands of zinc-blende BAs are shown in 
Fig. ||. Calculated band energies at high symmetry 
points are [gjvea-^f^Tablc VI.E-KJ As others have shown 
previously]i3't3'l 



acJ the band gap of BAs is indirect; the 
CBM is along the A line between the V and X points — 
at approximately 0.82(1, 0, 0)27r/a. An unusual feature of 
the band structure is the character of the CBM at T. In 
BAs, the CBM is the p-like T7 C state (ris c if spin-orbit in- 
teraction not included). Only the semiconductors silicon 
and BP share this feature. In most semiconductors, the 
lowest state at T is the singly degenerate s-like state. The 
origins of this feature of the BAs band structure can be 
understood in the context of the tight-binding model of 
HarrisonJl3 According to this model, the bonding (Fi5 V ) 



TABLE VI: Comparison of BAs band energies for 
LDA and GW.E3 Two LDA values are given: all- 
electron LAPW values from this, work and pseudopo- 
tential planewave (PP) values.Efl All the calculations 
used the experimental lattice constant of 4.777 A. The 
exchange-correlation potential used in the PP calcula- 
tions was not reported. For our LAPW calculations, 
the exchange-correlation of Perdew and Wangcil was used 
(See Sec. || for details.) 





LDA-PP (eV) 


LDA-LAPW (eV) 


GW (eV) 
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4.5 


4.57 
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r 8c (x2) 


3.3 
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r 7c 
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0.00 


0.00 


0.00 
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-0.22 
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-15.5 
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X% c 
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X-J C 


11.2 


11.13 
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Xdc 


1.38 
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1.86 
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X-J V 
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Xa v 


-11.3 
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-2.1 
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-8.82 
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Lev 


-12.6 


-12.62 


-13.6 



and antibonding (Pi5 C ) p states at P are given by 
E(T 15 ) ■- 




+ (4E„y (ii) 



where e c and e't are the p atomic orbital energies for the 
cation (c) and the anion (a) and the interatomic term 
E pp is proportional to l/%ond where dbond is the bond 
length. The bonding (Ti v ) and antibonding (Ti c ) s states 
at r are given by 



E(T 1 




(12) 



where e c s and e° are the s atomic orbital energies for the 
cation (c) and the anion (a) and the interatomic term E ss 
is again proportional to 1/^bond- A schematic diagram of 
the BAs energy levels in this model is shown in Fig. Ha. 
The reversal of the conduction band states in the 
BAs band structure relative to other III-V compounds 
is a result of two effects: (i) a small Tisv^TiSc bond- 
ing/ antibonding repulsion due to the small E pp , as de- 
scribed by Eq. ([ll]), and (ii) a large Tiv-Pic repul- 
sion of due to the small (ef — sf s ) 2 , as described by 
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FIG. 4: Schematic diagram of energy levels in the III— 
As series, (a) Schematic diagram of energy levels in zinc- 
blende BAs and their atomic origins. Because the B-p- 
As-p orbital energy difference is small relative to other 
III-V compounds and the B-s-As-s orbital energy dif- 
ference is large, the p-like r 15c conduction band state is 
below the s-like Ti c conduction band state. This ordering 
of the lower conduction band states is unusual for III-V 
compounds and is reminiscent of silicon, (b) Schematic 
diagram of the effects of p-d coupling on the VBM in 
AlAs, GaAs, and InAs. The VBM is driven up by the 
coupling of the 3d(4d) states in GaAs(InAs) and the As 
4p-like VBM, but in AlAs, this p-d repulsion drives the 
VBM down. 

Eq. (|l2|). As is evident from Fig. EL a small ri5 V -Ti5 C 
bonding/antibonding repulsion will lower the Ti5 C state 
and a large ri v -Ti c bonding/antibonding repulsion will 
raise the Tic. In BAs these two effects are enough to 
reverse the ordering of the Ti^ c and Tic states. 

The large ri v -Ti c bonding/antibonding difference is 
due to the large E ss due to the short bond length in 
BAs. Empirical tight binding has also shown that the 
r lv — Tic repulsion is much more|-p|ronounccd in BN, BP, 
and BAs than in other III-Vs.Ej The small Fi 5v -ri 5c 
bonding/antibonding repulsion is due to the small differ- 
ences in the p atomic orbitals of the boron and arsenic, as 
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FIG. 5: Partial and total densities of states for BAs 
[Eq. (g)] inside muffin tin spheres with radii 0.853 A for 
boron and 1.225 A for arsenic.EZI Note the strong mix- 
ing of B and As p states. Figure shows a quantitative 
measure of the p-p mixing in the each of the III-As com- 
pounds. 

shown in Fig. m. That is, the first term under the radical 
in Eq. ( [il| ) is small relative to other III- Vs. 

The strong hybridization of the B and As p states is 
visible in the partial densities of states plot [Eq. (g)] in 
Fig. g and band character plots in Fig. || [Eq. ([[)] where 
the width of the lines indicates the amount of the charac- 
ter indicated in the plot label. From these plots, one can 
see that the s band of arsenic between —12 and —16 eV is 
distinct and has very little mixing with other states. This 
is similar to the other three III-As systems, although in 
the case of BAs the arsenic s band is wider and some- 
what lower in energy. Two unique features of the BAs 
DOS compared to other III-V semiconductors is the pro- 
nounced boron p character at the VBM and anion p char- 
acter at the CBM. The stronger p-p hybridization in BAs 
relative to the rest of the III-As is due to the proximity 
of the orbital energies of B and As and the short bond 
length in the compound. 



C Trends in band gap and inter-valley energy 
differences along the III-As family 

Figure |7] shows the LDA-calculated vs. experimental 
band gaps for the III-As family, MAs (M = B, Al, Ga, 
In). For all but BAs, LDA energy band gaps are taken 
from the LAPW calculations reported in Ref. EJj. Exper- 
imental band gaps are from Ref. ES For BAs, for which 
no reliable experimental data exists, we used the GW 
band gap from Ref. p3. Because of the unusual ordering 





FIG. 6: Band characters of BAs [Eq. p)] plotted as a function of k. The thickness of the lines denotes a relative 
amount of a given character. The B and As p states are stongly mixed at both the VBM and the CBM whereas the 
As s band is distinct. There is some mixing of As p and B s states in the bottom of the p-p hybridized bands. 



of the conduction bands in BAs, the direct gap is not the 
Ti5v — * Tic gap but the Ti5 V — ► Ti5 C . The figures show 
the ri 5v — > Tic gaps (which are the direct gaps in AlAs, 
GaAs, and InAs). We find in Fig. [?] that the LDA vs. 
expt. Tisv — > Tic and the Fi5 V — + L c gaps lie almost 
on a straight line. The LDA error of these gaps in the 
III-As family is approximately constant, even for BAs. 
For the Ti$ v — ► Ti c gaps, the magnitude of the error is 
~1-1.5 eV. For the Fi 5v — > L c gaps, the magnitude of 
the error is ~ 0.7-1.1 eV. The experimental Ti$ v — > X c 
gaps only vary by ~0.5 eV through the series and cor- 
responding LDA gaps are nearly identical, resulting in a 
tight clustering of the Tisv — > X c points in Fig. m. 

In Fig. ||(a), we show the Fi 5v — > T lc , T 15v — ► L c , 
and ri5 V — ► X c gaps for the III-As family. We see: (i) 
the ordering of the conduction band states in BAs and 
AlAs is Ai c < Lie < Tic whereas (ii) the ordering is 
Tic < Lie < Ai c for GaAs and InAs, (iii) the X\ c and 



Lie states are closely spaced in AlAs and GaAs, (iv) BAs 
is "strongly indirect;" L\ c and Ti c are far above the CBM 
(which is close to Xi c ), and (v) a crossing of the lowest X 
and L conduction states occurs between AlAs and GaAs, 
but the Tisv — > Ti c gaps decrease so rapidly in the series 
that the Tic states lie below the L\ c in GaAs and InAs, 
and thus these two members of the family are direct band 
gap materials (whereas BAs and AlAs are indirect X gap 
materials). 

Figure H(b) shows the conduction band intervalley dif- 
ferences Tic — > X c and Ti c — > L\ c . Both differences show 
a smoothly decreasing trend from AlAs to GaAs to InAs, 
but an anomalously large decrease from BAs to AlAs. 
Again, this is a result of the enhanced s-s repulsion in 
BAs. Because the Ti c state has much more s character 
than either the lowest X and L conduction states (see 
Fig. o) , the effect of the s-s repulsion is greater for the 
Tic than for these other two states, driving it up relative 
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FIG. 7: Trends in the LDA vs. experimental band gaps of 
the III-As family. For the Ti$ v — > Ti c and the Ti5 V — > L c 
gaps, the errors are nearly constant. The LDA errors are 
~ 1-1.5 eV and -0.7-1.1 eV for the r 15v -» T lc and the 
ri5v — * -^c gaps, respectively. 

to them. This effect is more pronounced for the Ti c — ► X c 
difference because the lowest X conduction state has less 
s character than the lowest L conduction state. 
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FIG. 8: Trends in the inter-valley energy seperations of 
the III-As family. Note the crossing of the Ti5 V — > X c 
gaps and the ri 5v — > Ti c gaps between AlAs and GaAs. 
Thus, GaAs and InAs are direct gap materials whereas 
BAs and AlAs are indirect X gap materials. 



D Charge density and ionicity trends in the III-As 
family 

The valence charge density [upper four bands, d bands 
not included; Eq. (0)] in the (110) plane is shown for 
the four members of the III-As family in Fig. O. We see 
that the charge densities of AlAs, GaAs, and InAs are 
drawn towards the anion and exhibit a "single hump" in 
the bond charge. (The relative shifts in the bond charge, 
rfmax/4ond, are — 0.68 for AlAs and GaAs, and ~ 0.71 for 
InAs.) The relative shift of charge towards the anion (As) 
is due to the difference in electronegativity of the anion 
and the cation and reflects the partially ionic nature of 
the bond. In contrast, because of the almost entirely 
covalent nature of the bonding in the BAs system, in this 
system the charge density features a "double humpJi-band 
similar to the bond charges in diamond C and Si.c3e3 A 
similar analysis of the charge density .of BAs was reported 
earlier by Wentzcovitch and Cohen.Ej We see that BAs 
is unique in the III-As family in that the bonding is very 
covalent. This is evidenced by the exceptionally strong 
hybridization of p states from the anion and cation (see 



Figs. 0, p], and 0) and the significant bonding charge 
visible in the valence charge density which is similar to 
the bonding charge in diamond C and Si. 

It is interesting to note that, in contrast to AlAs, GaAs, 
and InAs where charge is drawn towards the anion, the 
nearly symmetric bond charge distribution shown in the 
valence charge plot of BAs is not centered exactly around 
the bond center but is actually drawn slightly toward the 
cation B. (The distance between the bond center and the 
boron position is about 45% of the total bond length.) 
This would indicate that whatever small ionic component 
exists in BAs results not from charge transfer towards 
the anion but from charge transfer towards the cation 
B. Further evidence of an ionic component in BAs where 
the boron atom acts as the anion is shown in the valence 
charge difference plots [Eq. (||)] in Fig. [L0| In AlAs, 
GaAs, and InAs, the maximum (marked as a cross) in 
Apvai occurs closer to the As anion than to the cation, 
consistent with the reported ionicities. However, in BAs, 
the maximum in Ap va \ occurs closer to the B atom. In 
addition, the positive region around the B atom in Ap va \, 
which does not occur for the other members of the III-As 
family, is further indication of a slight ionic component 
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FIG. 9: Valence charge densities p V ai [Eq. (Eh] for the III- 

As family, excluding the cation and anion d states. The 
contour plots use the same scale as the corresponding line 
plots. The contour spacing is 2 e/cell. The range of the 
line plots is from one tetrahedral interstitial site, through 
the bond, and to the neighboring tetrahedral interstitial 
site. Note the striking asymmetry of the bonding charge 
for the three systems AlAs, GaAs, and InAs, while the 
bond charge in BAs reveals a covalent, double-hump fea- 
ture similar to that found in diamond C, Si, or Ge. 



in BAs where B plays the role of the anion rather than 
the cation. This reversal of the cation and anion roles 
in BAs (and in BP-1 is predicted by the Phillips scale of 
electronegativitiescHI and has— been discussed previously 
by Wentzcovitch and Cohen.E] 



IV. BAND OFFSETS: BAs/GaAs/AlAs 

The offset AE V (AX/BY) between the valence band 
maxima of two semiconductor compounds AX and BY 
forming a heterostructure is one of the most important 
parameters in interfacial structures deciding both trans- 
port and quantum confinement. The results for the natu- 
ral (unstrained) band offsets [Eq. ra)] for BAs/GaAs and 
BAs/AlAs are plotted in Fig. [BJ The conduction band 



offsets shown in the figure are determined by adding the 
experimental gaps to the calculated valence band offsets. 
The natural bancLoffset for AlAs/GaAs has been com- 
puted previously^ and has been included here to con- 
trast the offset behavior of "typical" heterojunctions and 
those of heterojunctions with BAs. 

The ordering of the VBM of BAs in the III- As is some- 
what surprising. Because the VBM decreases systemat- 
ically in the order InAs -> GaAs — ► AlAs [Fig. [n|(b)], 
one might expect that the VBM of BAs would lie below 
AlAs, but in fact, the VBM of BAs is above AlAs, nearly 
as high as GaAs. The ordering of the VBMs for InAs, 
GaAs, and AlAs can be understood qualitatively using 
a tight-binding argument analogous to Eq. ( |TT| ) for ri 5c 
when a p-d term is added. The effects of this additional 
term are indicated schematically in Fig. 0b. The ordering 
of the VBMs is as follows: 

InAs relative to GaAs: Two effects account for the 
slightly higher VBM of InAs relative to that of GaAs.efl 
(i) InAs has a longer bond length than GaAs. This re- 
duces E pp of Eq. (^lj) and tends to drive the Fi 5v VBM 
up relative to GaAs. (ii) In tetrahedral symmetry, the 
cation d and anion p states share the same Tis symme- 
try and hence can interact through the potential matrix 
element ((/)i ation \V\(f> p amon ) = E dp . The mostly anion p- 
like VBM is repelled upwards by the Ga/In d stat es by . 
an amount E% p /(e% nion - £ c d aUon ). This p-d repulsionBEl 
is slightly stronger in InAs where the cation Ad states 
are shallower and less localized than the Ga 3d state, as 
shown schematically in Fig. [|b. 

GaAs relative to AlAs: The difference between the 
VBM of GaAs and AlAs can also be understood in terms 
of this p-d repulsion effect. In GaAs, the deep gallium 
d states couple to the arsenic p states, driving the VBM 
up, but in AlAs, there are no aluminum d states below 
the VBM. Instead, the interaction of the arsenic p states 
at the VBM and the unoccupied aluminum d states above 
the VBM drives the AlAs VBM down. Thus, p-d repul- 
sion has the opposite effect on the VBM in GaAs and 
AlAs (Fig. |b). 

AlAs relative to BAs: Two effects increase the VBM 
of BAs relative to the VBM of AlAs. (i) Because the 
unoccupied d states in BAs lie very high in energy rel- 
ative to the VBM, the p-d repulsion effect that drives 
the VBM down in AlAs is weaker in BAs. (ii) The more 
important effect is the unusual character of the VBM in 
BAs. In the rest of the III- As family, the character of the 
VBM is primarily anion p-like, but in BAs the bonding is 
much more covalent and the VBM comes from both the 
anion and the cation (Fig. |6|). Since the cation p levels 
are higher in energy than the As p levels (Fig. [y), any ad- 
mixture of cation p character pulls the VBM up. These 
two effects, lack of p-d repulsion and strong hybridiza- 
tion of the cation p character into the VBM, account for 
the high VBM of BAs relative to AlAs. 

GaAs relative to BAs: While the VBM of BAs and 
GaAs both lie above that of AlAs for the reasons given 
above, the VBM of BAs lies only slightly below that of 
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GaAs despite the much smaller lattice constant (~17% 
mismatch) and the lack of p-d repulsion. This is due to 
the unusual cation p hybridization into the VBM of BAs 
which raises the VBM by several tenths of an eV. 

Finally, we note that the transitivity among the band 
offsets of BAs/AlAs/GaAs is similar to other compounds 
where the lattice mismatches are not as large. By transi- 
tivity we mean that the band offset AE V (A/C) between 
two compounds A and C is well approximated as the 
sum of the band offsets between compounds A and B 
and between compounds B and C, i.e., 



We also modeled a 50%-50% random alloy using a 
32-atom-£pecial quasirandom structure (SQS16 in Ta- 
ble p).ta We found that the distribution of the bond 
lengths has_the expected bimodal form for a random bi- 
nary alloy,E3 and the B-As bonds are generally larger 
than the ideal B-As bond length while the Ga-As bond 
lengths are smaller than the ideal Ga-As bond length. 



AE v (GaAs/BAs) ^ AE v (BAs/AlAs) 



AE V (AlAs /GaAs). 
(13) 



In this case, AE v (GaAs/BAs) calculated directly yields 
0.19 eV but AE v (BAs/AlAs) + AE v (AlAs/GaAs) yields 
0.12 eV, a non-transitivity difference of 0.07 eV. 



V. BAs-GaAs ALLOYS 
A Bond lengths and bond angles in the alloys 

Interest in BAs-GaAs alloys centers around the hope 
that boron will modify the GaAs band gap similarly to 
nitrogen without the adverse effectEj of introducing lo- 
calized states that reduce carrier diffusion length. We 
know from the theory of III-V alloyscjEi that the opti- 
cal properties are decided by both the atomic relaxation 
and by charge transfer. Thus, we first study the bond 
relaxation around a boron substitutional impurity. 

The bond lengths and angles in the (110) plane of 
GaAs with ^3% boron substitution (BiGa 3 iAs32 super- 
cell) are shown in Fig. [lj. Bond lengths are shown as a 
percentage of the pure bulk GaAs bond length, except for 
the B-As bonds where the values indicate a percentage 
of the pure bulk BAs bond length. The numbers that lie 
between a triplet of atoms indicate the bond angle be- 
tween the three atoms as a percentage of the ideal bond 
angle in the zinc-blende structure of 109.47°. The rect- 
angle represents the supercell boundary. We note several 
observations: 

(i) Bond angles near boron increase by up to 5% from 
their ideal value. This helps accomodate the Ga-As and 
B-As alloy bonds (hereafter RguAs and RbAs, respec- 
tively) to keep the bond lengths close to their ideal val- 
ues, R BAs and i?g Qj 4 s in the pure binary compounds. 

(ii) The B-As bond decreases from the GaAs value 
towards the BAs value and ends up only 4% higher than 
the value for pure BAs. |— , 

(iii) The average bond length relaxation parameteiEa 



e = [R GaAs (x) - R BAs (x)]/[R a GaAs - R 



o i 
BAs\ 



(14) 



is 0.76. (e = when there is no relaxation and 1 when the 
relaxation is full.) Thus, assuming RbAs to be Vegard- 
like (as in VCA) overestimates the bond length by ~13%. 
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FIG. 11: Valence and conduction band offsets (in eV) for 
BAs/GaAs/AlAs and for the III-As family. The valence 
band offsets are calculated directly using Eq. (Q). The 
conduction band offsets are obtained by adding the ex- 
perimental band gap values to the valence band offsets. 
Part (a) shows the both the valence band and conduction 
band offsets for BAs/AlAs, AlAs/GaAs, and GaAs/BAs. 
Both the direct Fi 5c gap (solid line) and the X ic gap 
(dotted line) are shown. The height of the BAs VBM 
was determined by averaging the valence band offsets of 
BAs with respect to AlAs and GaAs. Part (b) shows the 
valence band offsets in the III-As family. Surprisingly, 
the VBM of BAs is found to be above that of AlAs and 
not very far below that of GaAs. 
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FIG. 10: Deformation density Ap val (p s v °^ d - pf^ m ) plots for the III-As family in the (110) plane. Contours marked 
with dashed lines correspond to negative values of Ap va i. The units are e/A 3 and the contour spacing is 0.02. The 
maxima are marked by crosses. To determine Ap va i, the overlapping charge densities of free atoms are subtracted 
from the self-consistent charge density of the solid. The atomic configurations for the free atoms are d w s 2 p 1 and s 2 p 3 
for the cations and anions, respectively. The d bands of the cations, but not of the anions, are included in the valence 
bands. 



B Bowing in the dilute alloy 

Both measurements and calculations indicate that, for 

isovalent semiconductor alloys, the deviation of the band 

gap, AE g (x), from the average band gap, E g (x), of 

the cnustitiiputs is often well described by a quadratic 
termJ^Bbiy 

AE g (x) =bx(x- 1). 



For typical semiconductor alloys, the bowing parame* 
tcr b is normally less than 1 eV and is independentcj 
of the concentration x. However, in GaAsi_ x N x al- 
loys, where the lattice mismatch is large and the bond 
strength differences are significant, theJpowing parame- 
ter is strongly composition dependentE3 and can be as 
large as 20 eV. Due to the large lattice mismatch of BAs 
and GaAs and the strong B-As bonds, one might ex- 
pect that B^Gai-^As alloys would also show such large 



and composition-dependent bowing. We calculated the 
ri5v —* Ti c gaps in the experimentally-relevant composi- 
tion range of 0-10%. The bowing was calculated using 64 
atom supercells, as described in Sec. ||. In this composi- 
tion range, the direct Ti^ v — > Ti c gap is the smallest gap. 
The bowing was calculated using the LDA Ti^ v — > Ti c 
gaps of BAs, GaAs, (4.64 and 0.28 eV, respectively) and 
the LDA gaps at two alloy compositions, 3% and 6%. At 
the 3% composition, the bowing is b — 3.4 eV and at 
6% the bowing is b = 3.6 eV. Thus, the band gap bow- 
ing is ~3.5 eV and relatively composition-independent in, 
this composition range. A recent experimental study,E-l 
estimated the bowing to be b = 1.6 ± 0.3 or 2.3 ± 0.3 us- 
ing theoretical estimates for the direct gap of cubic BAs 
of 3.56LJ and 4.23.E3 However, because of the LDA er- 
ror in the band gaps (see Fig. ^), these estimated band 
gaps are too small, and thus, the bowing is underesti- 
mated. Because the LDA errors for the band gaps are 
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FIG. 12: Bond lengths and angles in BGa3iAs32 in the 
(110) plane. Bond lengths are shown as a percentage of 
the ideal GaAs bond length, except for the B-As bonds 
where the values indicate a percentage of the ideal BAs 
bond length. The numbers that lie between a triplet of 
atoms indicate the bond angle between the three atoms 
as a percentage of the ideal bond angle in the zinc-blende 
structure of 109.47°. The rectangle represents the super- 
cell boundary. 



approximately constant for the III-As series, the bow- 
ing can be reliably calculated if all the points (Ea, Eb, 
and Ea 1 _ :c b :c ) are taken consistently from LDA calcu- 
lations (as we did) because only the relative values are 
important. For computing the band gap bowing from 
experimental alloy data, the GW estimate of the BAs 
direct band gap (5.5 eV),c3 which corrects LDA errors, 
is a more appropriate value to use to calculate the bow- 
ing. Using the experimental alloy data from Ref. 11^, the 
GW-estimated band gap for BAs (5.5 eV), and the ep. 
perimental value of 1.42 eV for GaAs, Geisz estimatesEll 
the revised bowing parameter to be ~3.5 eV in excellent 
agreement with the theoretical prediction. 

The bowing of B x Gai_ x As is much smaller than the 
bowing for GaAsi-^N^ alloys in the same range. Accord- 
ingly, we predict that addition of boron to GaAs will raise 
the gap. For example, the predicted gaps for 2%, 3%, 5% 
or 10% concentrations of boron, the band gaps would be 
1.43, 1.44, 1.46, and 1.51 eV, respectively. Thus, the ad- 
dition of boron to GaAs will not lead to a 1 eV gap alloy 
as seen with nitrogen addition. We also find that the 
bowing is composition independent. This implies that 
the effects of boron on the band gap in GaAs are much 
less pronounced than those of nitrogen. This is not com- 
pletely unexpected, however, as it has been previously 
observed that the bowing is smallepin mixed cation sys- 
tems than in mixed anion systems SB 

We have also calculated the band gap of the BGaAs2 
CuPt-ordered compound finding that it is 0.12 eV below 
the 50%-50% random alloy modeled via the SQS16. 



C Enthalpy of mixing and stability of superlattices 

Table M gives the formation enthalpies [Eq. (ft)] for 
(001) ordered superlattices of BAs/GaAs and the mixing 
enthalpies of B^Gai^^As random alloys. The formation 
energy of ordered superlattices of orientation G and pe- 
riod n is 



AH(n, G) = E{A n B n ; G) - ^(E A 



E B ) 



(15) 



can be seperated into two contributions S (i) The con- 
situent strain energy, i.e. the energy of the infinite pe- 
riod superlattice. This is the strain energy associated 
with deforming pure A and pure B into the in-plane lat- 
tice constant a of the superlattice, and relaxing them 
in the direction G. (ii) The interfacial energy, i.e., the 
difference between AH{n, G) and the consituent strain 
energy, is defined by 



AH(n,G) = 



2I(n,G) 



AE cs (a,G). 



(16) 



We calculated AEcs(a,G) for BAs and GaAs, deformed 
to the average lattice constant a along G — (001). This 
gave 148 meV/atom. From AH(n, G = 001) of Table 
and AEcs we calculated the interfacial energy / for 
n = 1, 2, and 4. We found I = 19, 23, and 26 meV, 
respectively. While a larger superlattice period n would 
be required to determine the converged value of /, we see 
that the BAs-GaAs interface is repulsive. This is why 
AH(n, G) of Table |V] decreases with n: larger n reduces 
the proportional effect of the interfacial repulsion. 

So far we dealt with (001) superlattices. The mono- 
layer (BAs) i/ (GaAs) i (111) is particularly interesting 
since in other III-V alloys (e.g., GaInP2) it appears as 
a spontaneously ordered stucture during alloying.E3 We 
find a very high AH(CuPt) of 108 eV/atom, suggesting 
thermodynamic instability. 

For the random alloys in Table |y|, all of the ex- 
cess enthalpies are positive. We compare the mix- 
ing enthaplies as obtained by the LDA (which includes 
both size-mismatched-induced strain effects and charge- 
tranfer "chemical" effects) and the valence force field 
(VFF) method (which includes only strain effects). The 
trends and magnitudes are similar, showing that the 
strain effects dominate AH. In Fig. [l| we compare the 
VFF-calculated mixing enthalpy of nitrogen in GaAs and 
contrast it with that of boron in GaAs, showing that 
the mixing enthalpy of GaAsi-^N^ alloys is much higher 
than for B^Gai-^As alloys. These results indicate that 
it may be easier to alloy boron with GaAs. This im- 
plies that the range of bulk B^Gai^As alloy composi- 
tions could be substantially larger than for the nitride 
alloy GaAsi-^Nj;. In epitaxial growth experiments, the 
alloy solubility can dramatically exceed that in bulk ex- 
periments for reasons explained in Ref. BQ. 
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FIG. 13: Calculated VFF mixing enthalpies for dilute 
GaAsi-zNj, alloys and dilute B^Gai-^As alloys. The 
mixing enthalpy is substantially smaller for boron alloys. 

D Quasi-localized electronic states 

We find that the incorporation of boron into the GaAs 
host material has little effect on the state at the VBM 
but that the lower conduction band states are strongly 
peturbed. The square of the wavefunctions for the VBM 
and the lowest two conduction band states for isolated 
boron in GaAs are shown in the left hand side of Fig. 
|l4| . The first column shows the states for an isolated 
boron atom in the center of a supercell of GaAs. The 
VBM is completely delocalized and looks very much like 
the VBM state of pure GaAs except for a region near 
boron impurity where the wavefunction is extended to- 
wards the boron atom. This distortion is precisely what 
one would expect based on the strong coupling of boron 
p and ar senic p sta tes at the VBM in BAs, as discussed 
in Sees. |III Bl and UlIDl. In contrast to the VBM, the 



CBM state shows a significant localization around the 
boron atom. The CBM shows long range delocalization, 
but the majority of the wavefunction is concentrated near 
the boron atom. For the second lowest conduction band 
state (CBM+1), the situation is similar except the wave- 
function is concentrated around a small number of gal- 
lium atoms as well. The dual character of the conduction 
band states (extended at long range but concentrated 
around the boron atoms in the short range) indicates 
that the states are resonant in the conduction band and 
are not localized states inside the gap. The conduction 
band states could be considered as boron-perturbed bulk 
GaAs conduction band states. 

The second column of Fig. [l4| shows the same states 
(VBM, CBM, CBM+1) for a fourth neighbor pair of 
boron atoms inside a supercell of GaAs. A statistically 
random distribution of boron atoms would result, among 
others, in pairs. In the case of GaAsi-^N^ alloys, the 
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FIG. 14: (color) Isosurfaces of the wavefunction squared 
for a single impurity and for a fourth nearest neighbor 
pair. Boron atoms are red, gallium black, and arsenic 
yellow. The isosurface value (blue) is 0.01 e/A 3 . In both 
cases, the VBM states are completely delocalized and are 
primarily As-derived. The lower conduction band states 
are strongly concentrated around the boron atoms. 

presence of nitrogen pair«-|Can result in localized impu- 
rity states inside the gap.Eil Calculations were performed 
for the five symmetrically-unique pair arrangements in 
the 64 atom supercell as described in Table [II . Qualita- 



tive features of wavefunction localization were the same 
for all of the pairs, and only the most representative case, 
the fourth nearest neighbor pair, is shown in the figure. 
The features of the pair states are similar to those of the 
isolated impurity discussed above — that is, the VBM is 
mainly an As-derived, delocalized state while the lower 
CB states are concentrated around the boron atoms. 
Again, there is some effect on the wavefunction very near 
the boron atom. The CBM state is concentrated around 
the boron atom in the middle of the cell and around the 
boron and gallium atoms on the cell edges. The next low- 
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est conduction band state shows similar features and the 
wavefunction is restricted primarily to the (200) plane of 
the cell. 

The short-range "localization" effects of boron incor- 
poration into GaAs appear to be similar to those seen in 
GaAsi_a;N x alloys, resulting mainly in "dual character" 
conduction band states that are still extended at long 
range but are localized around the boron atoms. How- 
ever, it appears that the perturbation of boron on the 
near-gap states of GaAs is gentler than that of nitrogen 
as none of the pairs cause states inside the gap. This 
is another instance (as with band gap bowing and band 
offsets) where boron is less "intrusive" than its first-row 
cousin, nitrogen. 



VI. CONCLUSIONS 

This paper explored the effects of alloying a conven- 
tional III-V compound, GaAs, with boron. Very little 
has been known, experimentally or theoretically about 
boride III-V semiconductor alloys, but boride semicon- 
ductor alloys are generating renewed interest now that 
epitaxial techniques have made it somewhat easier to fab- 
ricate them. In order to understand the fundamentals of 
this new class of possible alloys. We have studied zinc- 
blende BAs and its alloy with GaAs. 

For BAs we find: (i) The bonding in BAs is much 
more covalent than in the rest of the III-As family or 
other III- Vs. This is due to the similar electronegativi- 
ties of boron and arsenic, (ii) Consequently, BAs has a 
nearly symmetric, "double-hump" bond charge density, 
similar to silicon. The small asymmetry that does ex- 
ist in the bonding charge density actually favors boron. 
That is, charge is drawn toward the boron atom. In this 
sense, boron behave more like the anion than the cation 
in BAs. (iii) The band structure of BAs is more reminis- 
cent of silicon than most other III-V compounds. This is 
a consequence of (1) the small repulsion between bonding 
and antibonding p states and (2) the large repulsion of 
the arsenic and boron s states, (iv) The LDA errors in 
the band gaps are approximately constant in the III-As 
series BAs, AlAs, GaAs, and InAs. (v) The band off- 
sets of BAs/GaAs and BAs/AlAs are unexpected. The 
VBM of BAs is 0.39 eV higher than AlAs and only 0.19 
cV below GaAs. The main reason for this is the strong 
hybridization of both the cation and anion states at the 
VBM. 

For the B^Gai-^As alloys, we find that: (i) The bond 
angles are highly distorted near the boron sites. This ac- 
comodates the Ga-As and B-As bonds to stay close to 
their ideal bulk values. As a result, the bond length dis- 
tribution in a random alloy exhibits the typical bimodal 
behavior where each mode is near the individual ideal 
bulk values, (ii) The bowing parameter of B^Gai-^As 
for low concentrations of boron is 3.5 eV, much smaller 
than for GaAsi^Nj, alloys in the same composition 



range. Consequently, unlike GaN, alloying GaAs with 
small amounts (< 10%) of BAs increases the gap. (iii) 
The enthalpies of mixing indicate that the bulk solubility 
of boron in GaAs may be significantly higher than nitro- 
gen in GaAs, possibly increasing the composition range 
over which boride alloys may be formed, (iv) The wave- 
function of the VBM is completely delocalized and re- 
tains the character of pure GaAs. The lower conduction 
band states exhibit a "semi-localized" behavior — i.e. the 
states are strongly concentrated around the boron atoms, 
but far away they are extended states. 

In summary, the behavior of B^Gai^^As alloys is qual- 
itatively different from GaAsi_ x N x alloys. The perturba- 
tion of boron in GaAs is much "gentler" on most features 
of the electronic structure of GaAs than is nitrogen. If the 
same behavior is manifest for other boride III-V alloys, 
boron offers a new class of III-V alloys to be explored for 
novel behavior and possible device applications. Boron 
may also play an important role in alloys where it could 
be used as a relatively benign component which is added 
to lattice match the alloy to a given substrate. 
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